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S1-S3 counter charges in the voltage sensor module of a mammalian 
sodium channel regulate fast inactivation 



James R. Groome and Vern Winston 

Department of Biological Sciences, Idaho State University, Pocatello, ID 83209 



o 



The movement of positively charged S4 segments through the electric field drives the voltage-dependent gating 
of ion channels. Studies of prokaryotic sodium channels provide a mechanistic view of activation facilitated by 
electrostatic interactions of negatively charged residues in SI and S2 segments, with positive counterparts in the 
S4 segment. In mammalian sodium channels, S4 segments promote domain-specific functions that include activa- 
tion and several forms of inactivation. We tested the idea that S1-S3 countercharges regulate eukaryotic sodium 
channel functions, including fast inactivation. Using structural data provided by bacterial channels, we constructed 
homology models of the S1-S4 voltage sensor module (VSM) for each domain of the mammalian skeletal muscle 
sodium channel hNa v 1.4. These show that side chains of putative countercharges in hNa v 1.4 are oriented toward 
the positive charge complement of S4. We used mutagenesis to define the roles of conserved residues in the extra- 
cellular negative charge cluster (ENC), hydrophobic charge region (HCR), and intracellular negative charge clus- 
ter (INC). Activation was inhibited with charge-reversing VSM mutations in domains I— III. Charge reversal of ENC 
residues in domains III (E1051R, D1069K) and IV (E1373K, N1389K) destabilized fast inactivation by decreasing 
its probability, slowing entry, and accelerating recovery. Several INC mutations increased inactivation from closed 
states and slowed recovery. Our results extend the functional characterization of VSM countercharges to fast in- 
activation, and support the premise that these residues play a critical role in domain-specific gating transitions for 
a mammalian sodium channel. 
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INTRODUCTION 

Voltage-gated sodium channels initiate the action po- 
tential in excitable tissues (Armstrong and Hille, 1998; 
Catterall, 2012). These channels respond to membrane 
depolarization by activating, followed rapidly by an 
inactivating transition that limits the duration of the 
action potential (Hodgkin and Huxley, 1952). The gene 
sequence for the sodium channel reveals a pattern of 
regularly spaced arginine and lysine residues in con- 
served S4 transmembrane segments (Noda et al., 1984) . 
Channel gating is perturbed with mutations of these 
positively charged residues, suggesting that S4 segments 
act as voltage sensors to promote activation and initi- 
ate fast inactivation (Stiihmer et al., 1989). Additional 
experimental approaches have supported the premise 
that S4 segments carry the gating charge that translo- 
cates outward in response to membrane depolarization 
in sodium and other voltage-gated ion channels (Shao 
and Papazian, 1993, Perozo et al., 1994; Yang et al., 1996; 
Horn et al., 2000). Domain-specific roles for sodium 
channel S4 segments are suggested by experiments dem- 
onstrating that S4 translocation is coupled to fast inacti- 
vation (Chahine etal., 1994; Chen etal., 1996; Cha etal., 
1999; Sheets et al., 1999), to an intracellular activation 



gate (Muroi etal., 2010) , and to the selectivity filter (Capes 
etal., 2012). 

Electrostatic interactions of S4 positive charges with 
countercharges in surrounding segments SI, S2, and S3 
have been proposed as the stabilizing mechanism fa- 
cilitating S4 segment translocation. In this structural 
view, S1-S4 segments comprise a voltage sensor module 
(VSM) that surrounds a central, S5-S6 pore module 
(Chanda and Bezanilla, 2008). Experimental evidence 
for functional interactions within the VSM was first de- 
scribed for Shaker potassium channel (Papazian et al., 
1995) and subsequently extended to explain gating tran- 
sitions in Shaker and other potassium channels (Tiwari- 
Woodruff et al., 2000; Campos et al., 2007a; Zhang et al., 
2007; Wu et al., 2010; Pless et al., 2011). In sodium chan- 
nels, specific interactions of S4 residues with negative 
countercharges in SI and S2 segments have been identi- 
fied in the prokaryotic channel NachBac (DeCaen et al., 
2008, 2009, 2011). There, interactions within the VSM 
define the resting, closed state of the channel (Paldi and 
Gurevitz, 2010), with sequential ion pair interactions 
during depolarization that stabilize a series of activated 
states (Chakrapani et al., 2010; Yarov-Yarovoy et al., 2012) . 
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In this work, we sought to determine the functional 
roles for putative countercharges in the VSM of a 
mammalian sodium channel, hNayl.4. Eukaryotic volt- 
age-gated sodium channels are comprised of four ho- 
mologous domains with asymmetric S4 charge content 
(Trimmer et al., 1989) . We wished to identify the contri- 
butions of negative countercharges in S1-S3 segments 
of hNayl.4 toward domain-specific channel functions. 
To do this, we built sequence alignments and constructed 
homology models for each of four VSMs in hNayl.4, 
based on the structural data from bacterial sodium 
channels Na v Ab (Payandeh etal., 2011, 2012) andNa v Rh 
(Zhang et al., 2012). Our models predict that for each 
domain, negative countercharges in SI, S2, and S3 seg- 
ments are positioned appropriately to interact with pos- 
itive charges in the S4 segment. We then tested the role 
of VSM countercharges in hNayl.4 using site-directed 
mutagenesis and voltage clamp electrophysiology. We 
found that charge-reversing mutations in S1-S3 segments 
alter sodium channel activation and fast inactivation in a 
manner consistent with our hypothesis that negative 
countercharges are an important biophysical substrate of 
domain-specific voltage gating in the mammalian sodium 
channel. Portions of this work have appeared in abstract 
form (Groome, J.R., and V. Winston. 2012. Biophysical 
Society 56 th Annual Meeting. Abstr. 425b) . 

MATERIALS AND METHODS 

Alignments and homology modeling 

The crystal structure of the bacterial voltage gated sodium chan- 
nel Na v Ab (Payandeh et al., 2011) was downloaded from the 
Protein Database site (accession no. 3RVY; http://www.rcsb.org). 
Using DeepView (Swiss pdb viewer http:/ /spdbv.vital-it.ch/), the 
VSM of the A chain of 3RVY.pdb was saved as a separate file 
(3RVY_A_sensor.pdb), used as a template for subsequent model- 
ing. The amino acid sequence of the a subunit of human sodium 
channel protein type 4 (NP_000325) was downloaded from 
GenBank. Sequences of VSMs for domains I-IV were identified 
and saved as individual FASTA formatted files. 

Domain I-fV sequences were aligned to each other and to 3RVY_ 
A_sensor.pdb, using multiple sequence alignment in MUSCLE 
(http://www.drive5.com/muscle/). Each of the hNa v 1.4 VSM 
sequences was paired with the aligned amino acid sequence of 
3RVY_A_sensor.pdb, using the PIR format required by Modeller. 
Using the Modeller script model-single. py, a homology model 
structure was predicted for each hNayl.4 VSM. "Best" structures 
were identified using Modeller-derived GA341 scores and aligned 
with bacterial structures using the VMD multiseq plugin. 

Homology models of the VSMs in Na v Ab(2) (inactivated state 
of Na v Ab; Payandeh et al., 2012) and Na v Rh (NachBac orthologue; 
Zhang et al., 2012) were constructed using the same approach as for 
Na v Ab. Template structures were derived from 4EKW.pdb (acces- 
sion no. 4EKW; Na v Ab(2)) or 4DXW.pdb (accession no. 4DXW; 
Na v Rh) to construct homology models of the VSM for hNa v 1.4 in 
each domain. 

Molecular biology 

Mutations of the a subunit of human skeletal muscle sodium 
channel SCN4A were created with QuikChange XL II Site Di- 
rected Mutagenesis kits (Agilent Technologies) by incorporating 



mutant oligonucleotides into template SCN4A in vector pSP64T. 
Messenger RNA for a subunit was constructed from template 
linearized with EcoRl, using SP6 mMESSAGE m MACHINE 
in vitro transcription kits (Ambion). For all constructs, we co- 
expressed a subunit with subunit that associates with sodium 
channels and modulates fast inactivation (Yang et al., 1993). 
P subunit in vector pgH19 was linearized with Hindlll, and tran- 
scribed with T3 RNA polymerase. mRNA was injected into Xeno- 
pus laevis oocytes at a ratio of 3:1 Pi to a subunit in a volume of 
50 nL. Animals were maintained and surgical procedures were 
performed using protocols approved by the IACUC committee at 
Idaho State University. Oocytes were maintained in sterile media 
containing 96 mM NaCl, 2 mM KC1, 1 mM MgCl 2 , 4 mM CaCl 2 , 
5 mM Hepes, 5 mM sodium pyruvate, 3% horse serum, and 
100 mg/liter gentamicin, pH 7.4, with gentle agitation at 17.5°C, 
for 3-7 d before recordings. 

Electrophysiology 

Experiments were performed at 15°C, in the cut-open oocyte 
clamp configuration, with a CA-1B amplifier (Dagan Corpora- 
tion) . External solution consisted of 120 mM iVmefhyl i>glucamine, 
2 mM Ca(OH) 2 , and 5 mM Hepes, pH 7.4. Internal solution con- 
sisted of 120 mM jty-methyl D-glucamine, 2 mM EGTA, and 5 mM 
Hepes, pH 7.4. Before voltage clamp, the oocyte was permeabi- 
lized with 0.3% saponin to provide electrical access. Membrane 
potential was held at — 1 20 mV, with capacitance and leak subtrac- 
tion from a p/8 protocol, using Pulse 8.67 (HEKA). Sodium 
currents were analyzed with PulseFit 8.67 (HEKA) and Igor Pro 
6.02 (WaveMetrics). Statistical analyses were done with Instat 
2.0 (GraphPad Software) using Student's t tests of equilibrium 
parameters (Table 1), and of kinetic parameters at a given volt- 
age (Table 2) . 

We measured activation of hNa v 1.4 or mutant channels depo- 
larized with 20-ms test pulses to voltages ranging from —90 mV to 
60 mV. Current amplitude in response to each test pulse was nor- 
malized to maximum, and the resulting curve was fit with a 
Boltzmann distribution according to the equation I/Imax = 1/(1 + 
exp[zF(V 0 .5 - V M )/RT]), where I is the response to a test pulse, 
Imax is the maximum response, z is the slope factor, F is Faraday's 
constant, V 0 . 5 is the midpoint voltage, V M is the command voltage 
for a given test pulse, R is the universal gas constant, and T is tem- 
perature in degrees Kelvin. From this equation, we derived pa- 
rameters of midpoint (V0.5) and slope factor of activation, and 
calculated the free energy change AG as 4zFVo.s according to a 
four state model as described in Chowdhury and Chanda (2012). 
We measured activation kinetics as 10-90% rise time to peak cur- 
rent over a voltage range from — 20 mV to 60 mV. 

We determined steady-state fast inactivation by precondition- 
ing oocyte membrane potential for 300 ms to voltages ranging 
from —120 mV to 25 mV and then delivering depolarizing test 
pulses to 30 mV. Peak current amplitude in response to each test 
pulse was normalized against the maximum, and the resulting 
h„ curve was fit with a Boltzmann distribution according to the 
following equation: I/Imax = 1 ~~ 1/(1 + exp[zF(h 05 - V M )/RT]), 
where V M is the conditioning voltage, to derive parameters of mid- 
point (ho.5) and slope factor for inactivation. Kinetics of entry into 
fast inactivation were obtained from exponential fits to decays 
of sodium current over the voltage range of —90 mV to 60 mV ac- 
cording to the equation I, = Ko + Ki exp( — t/taun), where I is so- 
dium current at time t, Ko is the asymptote, Ki is initial amplitude 
of sodium current, and tau H is time constant of fast inactivation. 

Parameters for recovery of channels from fast inactivation were 
determined using a double pulse protocol (Fig. SI). Channels 
were inactivated with a 30-mV depolarization for 30 or 50 ms, fol- 
lowed by command hyperpolarization for durations up to 53 ms, 
at voltages ranging from — 1 20 mV to — 70 mV, and test pulses to 
30 mV. The normalized recovery curve was fit with the previous 
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TABLE 1 

Equilibrium parameters for wild-type and mutant sodium channels 



Mutation 


Location 


v 

V 0.5 


Slope factor 


AAG° 


u 

"0.5 


Slope factor 






mV 




heal/ mol 


mV 




hNa v 1.4 




18.3 ± 1.2 


1.39 ±0.03 




-65.8 ± 1.1 


4.24 ±0.12 


ENC 














S151K 


Domain I SI 


23.2 ± 0.8 b 


1.37 ±0.06 


0.6 


-66.5 ± 0.7 


4.67 + 0.06 b 


D152R 


Domain I SI 


23.2 ± 0.9 b 


1.37 ±0.03 


0.6 


-64.1 ±0.6 


4.11 ± 0.25 


E161R 


Domain I S2 


23.8 ± 1.0 C 


1.57 ± 0.03 d 


1.1 


-60.5 ± 0.6 C 


3.82 + 0.08 b 


E598K 


Domain II SI 


36.3 ± 0.9 d 


1.59±0.05 b 


3.0 


-59.8 ± 1.4 b 


2.90 ± 0.14 d 


N614K 


Domain II S2 


30.3 ± 0.8 d 


1.52±0.03 b 


1.9 


-61.0±0.6 d 


3.84 ± 0.20 


E1051R 


Domain III SI 


26.0 ± 0.7" 


1.60 ±0.02" 


1.5 


-48.4 ± 1.2 d 


3.16±0.16 d 


E1051D 


Domain III SI 


25.7 ±1.3" 


1.39 ±0.03 


1.0 


-62.1 ± 1.2* 


.3.98 ± 0.15 


D1069K 


Domain III S2 


30.0 ± 0.8 d 


1.61 ±0.04 c 


2.1 


-49.3 ± 1.9 d 


.3.56 ± 0.1.3 C 


D1069E 


Domain III S2 


19.2 ± 1.2 


1.43 ±0.02 


0.2 


-64.4 ±0.7 


2.88±0.11 d 


E1373K 


Domain IV SI 


15.2 ± 1.0 


1.73 ±0.03" 


0.1 


-50.3 ± 0.7 d 


2.47 ± 0.07 d 


E1373N 


Domain IV SI 


18.1 + 1.2 


1.60±0.02 d 


0.3 


-49.9 ± 0.9 d 


2.41 ± 0.06 d 


E1373D 


Domain IV SI 


20.4 ± 0.9* 


1.39 ±0.02 


0.3 


-65.8 ± 0.7 


2.81 ± 0.07 d 


N1389K 


Domain IV S2 


18.6 ± 1.1 


1.65 ± 0.03 d 


0.5 


-48.7±0.9 d 


2.58 ± 0.12 d 


N1389E 


Domain IV S2 


21.8 + 1.0* 


1.31 ±0.01" 


0.3 


-71.5 + 1.0° 


4.59 ±0.11* 


N1389D 


Domain IV S2 


19.2 ±0.7 


1.33 ±0.02 


0.0 


-75.4±0.6 d 


2.66 ± 0.03 d 


HCR 














N144K 


Domain I SI 


35.4 ± 1.2" 


1.58 ±0.07* 


2.8 


-61.7±0.9 b 


3.69 ± 0.70 b 


N144R 


Domain I SI 


34.2 ± 1.2" 


1.40 ±0.03 


2.1 


-62.3 + 0.5 b 


4.06 ± 0.09 


N144D 


Domain I SI 


20.2 ±0.7 


1.19±0.01 d 


-0.1 


-72.0±0.5 d 


5.13±0.09 d 


N144E 


Domain I SI 


20.1 ±0.9 


1.29±0.02 b 


0.1 


-67.7 ± 0.3 


5.06 ± 0.09 d 


N591K 


Domain II SI 


34.3 ± 1.7" 


1.30 ±0.03* 


1.8 


-58.5 ± 1.0 d 


2.89 ± 0.1.3 d 


N591R 


Domain II SI 


34.4 ± l.l d 


1.38 ±0.03 


2.0 


-61.0 + 0.4" 


5.93 ± 0.15 d 


N591D 


Domain II SI 


32.0 ± l.l d 


1.50 ±0.04* 


2.1 


-61.2±0.4 C 


4.54 ± 0.07* 


N591E 


Domain II SI 


18.8 ± 1.3 


1.42 ±0.02 


0.1 


-63.8 ±0.7 


4.52 ±0.11 


S1044K 


Domain III SI 


23.6 ± 0.8 C 


1.43 ±0.03 


0.8 


-57.9 + 0.6 d 


4.00 ± 0.1.3 


S1044R 


Domain III SI 


30.4 ± 0.7 d 


1.42 ±0.02 


1.6 


-49.8 + 0.6 d 


4.22 ±0.10 


S1044N 


Domain III SI 


22.8 ± 0.6 b 


1.39 ±0.02 


0.6 


-60.6±0.6 C 


4.55 ± 0.07* 


S1044D 


Domain III SI 


21.4 ± 1.2 


1.32 ±0.02* 


0.2 


-70.9±0.5 C 


.3.99 ± 0.09 


N1366K 


Domain IV SI 


19.3 ± 1.6 


1.41 ±0.02 


0.2 


-64.2 ± 1.2 


.3.88 ± 0.08* 


N1366R 


Domain IV SI 


21.5 ± 1.4 


1.51 ±0.11 c 


0.7 


-58.3 ± 0.8 d 


2.64 ± 0.04 d 


N1366D 


Domain IV SI 


26.3 ± 1.7 C 


1.25±0.04 b 


0.7 


-85.8 ± 0.6 d 


3.37 ± 0.09 d 


N1366E 


Domain IV SI 


21.8 ± 1.0* 


1.41 ±0.02 


0.5 


-71.2 ± 1.3 b 


.3.44 + 0.08 d 


INC 














E171R 


Domain I S2 


36.7±2.1 d 


1.24±0.05 b 


1.9 


-75.7 ± 1.7 d 


3.79 ± 0.25 


K175E 


Domain I S2 


22.1 ± 1.2" 


1.37 ±0.02* 


0.5 


-65.0 ± 0.9 


4.24 ±0.14 


N194K 


Domain I S3 


19.1 ± 1.0 


1.33 ±0.02 


0.0 


-68.7 + 0.8* 


2.88 ± 0.04 d 


D197R 


Domain I S3 


35.1 ± 1.5 d 


1.28±0.04 b 


1.8 


-76.1 + l.l d 


4.44 ± 0.24 


E624K 


Domain II S2 


26.5 ± 1.2 d 


1.44 ±0.03 


1.2 


-63.7 ±0.7 


5.13±0.14 d 


K628E 


Domain II S2 


18.3 ± 0.9 


1.34±0.02 c 


-0.1 


-68.9 + 0.9* 


.3.95 ±0.11 


N643K 


Domain II S3 


17.7 ±0.4 


1.41 ±0.02 


0.0 


-61.2±0.8 b 


4.51 + 0.08 


D646K 


Domain II S3 


27.6 ± 1.2 d 


1.41 ±0.03 


1.2 


-64.2 ±0.6 


3.95 ± 0.13 


E1079R 


Domain III S2 


21.6 ±0.9* 


1.57±0.05 b 


0.8 


-63.4 ± 1.1 


3.86 ±0.16 


K1083E 


Domain III S2 


20.7 ± 1.0* 


1.39 ±0.02 


0.3 


-68.9±0.7 b 


3.71 ±0.11* 


C1098K 


Domain III S3 


23.1 ± 1.0 b 


1.27±0.02 c 


0.4 


-75.7±0.6 d 


3.03±0.11 d 


D1101V 


Domain III S3 


20.2 ± 1.3 


1.44 ±0.02 


0.3 


-63.4 ± 1.0 


4.54 ± 0.10 


E1399K 


Domain IV S2 


20.2 ± 0.8 


1.47 ±0.02* 


0.4 


-61.0±0.9 b 


.3.29 + 0.08 d 


K1403E 


Domain IV S2 


20.7 ± 1.6 


1.50 ±0.04 


0.5 


-63.9 + 1.0 


3.03±0.11 d 


N1417K 


Domain IV S3 


24.2 ± 0.9 C 


1.38 ±0.01 


0.7 


-69.4 ± 0.6 b 


3.77 ± 0.06 b 


D1420K 


Domain IV S3 


22.4 ± 1.2* 


1.57±0.04 c 


0.9 


-70.1 ±0.8 b 


2.83 ± 0.09 d 



Mean values ± SEM from Boltzmann fits of I/V relations in 10-24 experiments. The free energy change in wild-type and mutant channels is calculated 

from 4zFV 0 .5. The change in free energy of activation produced by mutation is given as AAG°. 

*P < 0.05. 

b P<0.01. 

C P< 0.001. 

d P< 0.0001. 
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TABLE 2 





Kinetic parameters for activation ( 10—90% time to peak current, 20 mV), time constants of entry (20 


mV), and recovery from fast 






inactivation (—100 m 


V), and delay in onset to recovery (— 100 mV) 




Mutation 


T f\ f 1 -1 1 1 ^ Tl 


± Ullt: Wj IJCrtJs. cLV- LI V LlLlUll 


Fast inactivation tau 


Recovery tau 


Recovery delay 






ms 


HIS 


HIS 


WIS 


hNa v 1.4 




1.05 ±0.03 


1.26 ±0.03 


7.15 ±0.54 


1.04 ± 0.05 


ENC 












S151K 


Domain I SI 


1.33 ± 0.03 d 


1.23 ±0.05 


11.3±0.61 d 


2.18±0.16 d 


D152R 


Domain I SI 


1.04 ±0.04 


1.15 ±0.06 


8.09 ± 0.48 


1.42 ±0.19* 


E161R 


Domain I S2 


1.18±0.03 b 


1.31 ±0.07 


5.34±0.25 b 


0.95 ± 0.06 


E598K 


Domain II SI 


1.09 ±0.07 


2.13±0.27 b 


5.17 ±0.48* 


1.09 ±0.10 


N614K 


Domain II S2 


1.15 ±0.06 


1.65±0.10 b 


5.50 ± 0.38" 


0.86 ± 0.05* 


E1051R 


Domain III SI 


1.55 ±0.14" 


3.25 ± 0.55 b 


3.61 ± 0.24 d 


0.79 ± 0.07 b 


E1051D 


Domain III SI 


1.05 ±0.04 


1.23 ±0.06 


5.51 ± 0.30* 


1.17 + 0.10 


D1069K 


Domain III S2 


1.18 ±0.09 


2.92 ± 0.46 b 


3.42 ± 0.16 d 


0.84 ± 0.05 b 


D1069E 


Domain III S2 


1.18±0.05 a 


1.94±0.16 b 


6.31 ± 0.36 


1.15 ± 0.09 


E1373K 


Domain IV SI 


1.65 ±0.07" 


5.51 ± 0.31 d 


3.42 ± 0.39 d 


0.44 ± 0.07 d 


E1373N 


Domain IV SI 


1.54±0.07 d 


4.94±0.18 d 


3.70 ± 0.23 d 


0.46 ± 0.02 d 


E1373D 


Domain IV SI 


1.24 ± 0.03 d 


2.74±0.22 d 


4.85 ± 0.22 c 


0.81 ± 0.05 b 


N1389K 


Domain IV S2 


1.85±0.12 c 


4.14±0.29 d 


4.31 ±0.42 c 


0.42 ± 0.02 d 


N1389E 


Domain IV S2 


1.26±0.04 d 


1.79±0.07 d 


8.36 ±0.55 


1.25 ± 0.07" 


N1389D 


Domain IV S2 


1.21 ±0.02 d 


1.52±0.04 d 


10.66 ±0.39 d 


1.21 ± 0.10 


HCR 












N144K 


Domain I SI 


1.65 ± 0.1 l d 


2.51 ± 0.30 c 


6.19 ±0.37 


1.46±0.48 c 


N144R 


Domain I SI 


1.54 ±0.02" 


1.70±0.08 d 


7.63 ± 0.30 


0.94 ± 0.25 


N144D 


Domain I SI 


1.02 ±0.01 


1.11 ±0.04 c 


9.60 ± 1.1* 


1.39±0.11 c 


N144E 


Domain I SI 


1.05 ±0.02 


1.15 ±0.04* 


7.19 ±0.60 


1.09 ±0.05 


N591K 


Domain II SI 


1.56±0.09 d 


2.62±0.21 d 


5.68 ± 0.31* 


0.97 ± 0.06 


N591R 


Domain II SI 


1.43±0.02 d 


1.33 ±0.04 


5.51 ± 0.28* 


0.88 ± 0.07 


N591D 


Domain II SI 


1.25 ± 0.03 d 


1.09±0.03 c 


6.90 ± 0.43 


1.00 ±0.06 


N591E 


Domain II SI 


1.15 ±0.04" 


1.11 ±0.03 b 


7.90 ± 0.59 


1.53±0.12 c 


S1044K 


Domain III SI 


1.20±0.03 b 


1.63 ±0.17* 


6.46 ± 0.41 


1.14 ± 0.04 


S1044R 


Domain III SI 


1.30 ± 0.03 d 


1.55 ± 0.03 d 


4.49 ± 0.25 c 


0.67 ± 0.03 d 


S1044D 


Domain III SI 


1.23±0.04 c 


1.84±0.1 d 


9.94 ± 0.33 d 


1.96 ± 0.13 d 


S1044N 


Domain III SI 


1.10 ±0.03 


1.27 ±0.04 


5.82 ± 0.33* 


0.95 ± 0.07 


N1366K 


Domain IV SI 


1.17±0.03 b 


1.83±0.10 d 


5.89 ± 0.38 


0.94 ± 0.05 


N1366R 


Domain IV SI 


1.91 ±0.07 d 


4.09±0.16 d 


14.62 ± 0.58 d 


0.40 ± 0.07 d 


N1366D 


Domain IV SI 


1.09 ±0.03 


1.60 ± 0.03 d 


13.02 ± 0.88 d 


2.99 ± 0.13 d 


N1366E 


Domain IV SI 


1.37 ± 0.03 d 


2.05 ± 0.08 d 


6.45 ± 0.35 


1.37±0.06 c 


INC 












E171R 


Domain I S2 


1.09 ±0.04 


1.41 ±0.06* 


9.59 ± 0.48 b 


2.20±0.21 d 


K175E 


Domain I S2 


1.13 ±0.02" 


1.45 ±0.09 


6.21 ± 0.36 


1.05 ± 0.06 


N194K 


Domain I S3 


1.52 ±0.04" 


4.64 ± 0.15 d 


6.18 ±0.30 


1.13 ±0.04 


D197R 


Domain I S3 


1.02 ±0.04 


1.56±0.08 b 


16.1 ± 1.43 d 


2.80 ± 0.37 d 


E624K 


Domain II S2 


1.26 ±0.04" 


1.57±0.09 b 


6.11 ±0.28 


0.91 ± 0.08 


K628E 


Domain II S2 


1.12 ±0.03 


1.45 ±0.07* 


6.72 ± 0.15 


1.09 ±0.04 


N643K 


Domain II S3 


1.24±0.04 c 


1.49±0.05 c 


6.04 ± 0.37 


0.96 ±0.11 


D646K 


Domain II S3 


1.50±0.04 d 


2.09 ± 0.32* 


7.13 ±0.30 


1.05 ± 0.09 


E1079R 


Domain III S2 


1.54 ± 0.05 d 


2.24±0.18 d 


7.13 ±0.36 


1.07 ±0.09 


K1083E 


Domain III S2 


1.16 ±0.04* 


1.31 ±0.05 


8.22 ± 0.33 


2.18 ± 0.25 d 


C1098K 


Domain III S3 


1.11 ±0.02 


1.36 ±0.04* 


10.4±0.60 c 


1.61 ± 0.11 c 


D1101V 


Domain III S3 


1.09 ±0.02 


1.23 ±0.07 


7.63 ± 0.38 


1.26±0.05 b 


E1399K 


Domain IV S2 


1.57 ±0.04" 


3.25 ± 0.08 d 


7.75 ± 0.91 


1.21 ± 0.12 


K1403E 


Domain IV S2 


1.19 ±0.06 


2.47 ± 0.30 b 


5.44 ± 0.22 


1.10 ±0.06 


N1417K 


Domain IV S3 


1.23 ± 0.03 d 


1.57±0.06 d 


7.50 ± 0.32 


1.25 ±0.08* 


D1420K 


Domain IV S3 


1.45 ± 0.05" 


3.52 ± 0.48 c 


7.32 ±0.31 


0.80 ±0.11" 



Mean values ± SEM from 10-24 experiments. 

*P < 0.05. 

b P<0.01. 

C P< 0.001. 

d P< 0.0001. 
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Figure 1. Homology modeling of sodium channel VSM. (A) Sequence alignments of VSM from bacterial channels Na v Rh and Na v Ab, 
and domains I-IV of mammalian sodium channel hNa v 1.4. (B) Homology models of each domain of hNa v 1.4 based on the Na v Ab tem- 
plate, showing relative positions of putative countercharges in SI, S2, and S3 segments, and positively charged residues through R4 in 
S4 segments. (C) View of homology models in B with S4 residues R1-R3 aligned in each, to show orientation of countercharges to the 
S4 segment. 



equation to obtain tau H for recovery. Delay in onset to recovery 
from fast inactivation was taken as time at which normalized re- 
covery current became positive (Kuo and Bean, 1994). 

Online supplemental material 

Electrophysiology. An illustration of methodology for calculating 
recovery time constant and its delay from the double pulse proto- 
col is provided in Fig. SI. Figs. S4-S6 provide the voltage depen- 
dencies of recovery delay for mutations tested in this work. 

Modeling. Template and domain specific VSM models are 
available as text files: (1) NaVAb: 3RVY_A„sensor, 3RVY_NaV_DI, 
3RVY_NaV_DII, 3RVY_NaV_DIII, 3RW_NaV_DIV, (2) NaVAb(2) 
4EKW_sensor, 4EKW_NaV_DI, 4EKW_NaV_DII, 4EKW_NaV_ 
Dili, 4EKW_NaV_DIV, and (3) NaVRh: 4DKW_sensor, 4DKW_ 
NaV_DI, 4DKW_NaV_DII, 4DKW_NaV_DIII, 4DKW_NaV_DIV. 

We include homology models of the VSM for each hNayl. 4 do- 
main, based on structural data provided by NayAb(2) and Na v Rh 
as Fig. S2. We also provide homology models based on Na v Ab that 
include residues not investigated in this study but relevant to this 
work, in Fig. S3. 

Online supplemental material is available at http://www.jgp 
.org/cgi/content/full/jgp.201210935/DCl. 



RESULTS 

Identification of putative countercharges in hNa v 1 .4 
Alignments of the VSM from domains I-IV of hNayl .4 
to bacterial sodium channels NayAb (Payandeh et al., 

2011) and NachBac orthologue NavRh (Zhang et al., 

2012) are shown in Fig. 1 A. Conserved sequences rep- 
resenting putative countercharges in SI, S2, and S3 
segments are outlined in boxes. These alignments iden- 
tified amino acid residues in hNayl -4 homologous to 
bacterial sodium channel VSM residues described in 
crystal structure, characterized for ion-pair interactions, 
or suggested by structural modeling but not yet func- 
tionally characterized. 

From the bacterial NayAb template, homology mod- 
els of the VSM for each domain of hNayl -4 were cre- 
ated. Fig. 1 B shows the position of each of the conserved 
residues in the mammalian VSM, with respect to posi- 
tively charged residues in S4. With the S4 segment 
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aligned along the R1-R3 backbone (Fig. 1 C) , side chains 
for mammalian putative countercharges in domains 
I-IV are each oriented similarly, facing toward the S4 
segment. Models created from the Na\Ab(2) and Na v Rh 
templates are shown in Fig. S2. The progressively ac- 
tivated state of the channels in these three models is 
observed with the position of R4 beneath (NayAb), tra- 
versing (NayAb (2) ) , and above (NavRh) the gating pore. 
The orientation of putative countercharges in the VSM 
was similar in each model, with the exception of the 
outer SI residue in the extracellular negative charge 
cluster (ENC) , which was observed in a nonstructured 
region in the NayRh model. 

The sequence alignments and resulting homology 
models for each domain suggested to us that conserved 
residues in S1-S3 might influence voltage-dependent 
functions of hNayl.4. Therefore, we characterized each 
of these residues with mutagenesis, using voltage clamp 
experiments. We compared steady-state and kinetic pa- 
rameters between wild-type hNayl.4 channels and mu- 
tations in the ENC, intracellular negative charge cluster 
(INC) , and a conserved charge just above the putative 
gating pore, or hydrophobic charge region (HCR). The 



position of that HCR SI residue relative to the gating 
pore is illustrated in Fig. S3. Also shown are positions 
of isoleucine, valine, and aromatic residues in each 
domain of hNayl.4, homologous to those proposed in 
NayAb to comprise the hydrophobic constriction site 
(HCS; Payandeh et al., 201 1; Yarov-Yarovoy et al., 2012) . 

ENC mutations 

In the bacterial sodium channel NachBac, negatively 
charged residues in the ENC pair sequentially with posi- 
tively charged residues in S4 during activation (DeCaen 
et al., 2008, 2009, 2011). Our models predicted ho- 
mologous ENC countercharges in each domain of the 
mammalian channel hNayl.4, with some question as to 
the roles of S151, D152, or both as countercharges in 
domain I. We compared the effects on channel function 
for charge-reversing mutations of the outer residue 
at SI (including S151Kand D152R) in domain I, E598K 
in domain II, E1051R in domain III, and E1373K in 
domain IV, and at S2 (E161R in domain I, N614K in 
domain II, D1069K in domain III, and N1389K in do- 
main IV) . For ENC residues in domains III and TV, we 
also tested the effects of charge-conserving mutations. 



1.0 



0.8- 



0.6 



~ 0.4 



0.2- 



0.0 



hNa„1.4 




♦ hNa„1.4 
O E1051R 

• E1051D 
D1069K 

■ D1069E 



* 

i 

t m 

.? 



I' 



-80 -60 



-40 -20 0 20 
Voltage (mV) 



40 60 



D 



♦ hNa v 1.4 
O S151K 
M D152R 
□ E161R 
A E598K 
N614K 



S W 

A. a 
♦ * 



*••« 



##1' 



-80 -60 -40 -20 0 20 40 60 
Voltage (mV) 



1.0 



0.6 



0.4 



0.2 



hNa„1.4 

E1373K 
E1373D 
E1373N 
N1389K 
N1389D 
N1339E 



0 ^ 0 ^ ^ t 



i 



-60 -40 -20 0 20 40 60 
Voltage (mV) 



Figure 2. Conductance in hNa v 1.4 and ENC mutations. (A) Traces for wild type and mutant channels in response to depolarizing com- 
mands to voltages ranging from —90 mV to +60 mV. I/V relations are shown for ENC mutations in domain I and domain II (B), 
domain III (C), and domain IV (D). Values represent mean ± SEM (error bars) from 10-22 experiments. 
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Mutations of ENC residues in domains I— III inhibited 
activation (Fig. 2 and Tables 1 and 2). In domain I, 
activation midpoint (V0.5) was depolarized by 5 mV with 
mutation at SI (S151K, D152R) or at S2 (E161R). Acti- 
vation was slowed by S151Kand E161R, and E161R also 
increased slope factor. Larger depolarizing shifts in I/V 
relations were observed for domain II ENC mutations 
E598K (SI; 18 mV) and N614K (S2; 12 mV). In domain 
III, significant depolarizing shifts of V0.5 were produced 
byE1051R (SI, 8mV) and E1051D (7 mV). E1051R in- 
creased slope factor and slowed rise time to peak ac- 
tivation. D1069K (S2) significantly depolarized V0.5 by 
12 mV and increased slope factor, whereas D1069E 
produced activation midpoint and slope factor similar 
to hNayl.4. For domain IV, charge reversal or substitu- 
tion did not affect the probability of channel opening, 
as I/V relations for El 373K/N/D (SI) orN1389K/D/E 
(S2) were similar to hNayl-4. However, each of these 
mutations significantly slowed rise time to peak activa- 
tion. Slope factor was increased by mutations E1373K, 
E1373N, and N1389Kand decreased by N1389E. 

Most ENC charge-reversing mutations produced a de- 
polarizing shift in the midpoint (ho.5) of the steady-state 
fast inactivation (h„) curve (Fig. 3 and Table 1). Mu- 
tations E161R (domain I S2), E598K (domain II SI), 
and N614K (domain II S2) depolarized midpoint by 
4—6 mV, and E598K also decreased slope factor (Fig. 3 A) . 
Charge-reversing mutations at SI or S2 in the ENC of 
domain III (E1051RandD1069K) or domain IV (El 373K 
and N1389K) produced the largest effects on steady- 
state inactivation, with 15-17-mV depolarizing shifts of 



the h« curve and significantly decreased slope factor 
(Fig. 3, B-D). Charge-substituting mutations E1051D 
and D1069E produced h 0 .5 values similar to hNayl.4. 
E1051D, but not D1069E, produced a slope factor simi- 
lar to the wild-type channel. At domain IV SI, aspartate 
substitution at El 373 produced h 0 .5 similar to hNayl.4, 
whereas asparagine substitution produced an effect simi- 
lar to E1373K, and both mutations significantly decreased 
slope factor. At domain IV S2, glutamate or aspartate 
substitution at N1389 produced a significant hyperpo- 
larizing shift of the h x curve. Slope factor for N1389E, 
but not N1389D, was similar to hNayl.4. These findings 
indicate that the charge content of domain III and do- 
main IV ENC residues is an important determinant of 
fast inactivation. 

Fast inactivation was slowed by reversal of charge in 
the ENC for domains II-IV (Fig. 4 and Table 2). In do- 
main II, E598K (SI) and N614K (S2) each slowed the 
entry of channels into fast inactivation, at voltages more 
positive than -20 mV (Fig. 4 B). In domain III, E1051R 
(SI) and D1069K (S2) slowed fast inactivation over this 
voltage range, and decreased its voltage dependence 
(Fig. 4, C and D). E1051D fully restored, and D1069E 
partially restored, wild-type fast inactivation kinetics. 
In domain IV, charge reversal at SI (E1373K) or at S2 
(N1389K) dramatically slowed fast inactivation at all 
voltages tested (Fig. 4, E and F). E1373N slowed fast in- 
activation similar to E1373K, whereas El 373D produced 
kinetics intermediate between E1373K and hNayl.4. At 
domain IV S2, N1389D/E slowed fast inactivation to a 
lesser extent than N1389K. 
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Figure 3. Steady-state fast inactivation 
in hNayl.4 and for ENC mutations. 
Channels were conditioned to volt- 
ages ranging from — 120 mV to +25 mV 
for 300 ms before test pulses to 30 mV. 
I/V relations (h„ curves) are shown for 
ENC mutations in domain I and do- 
main II (A), domain III (B), domain IV 
SI (C), and domain IV S2 (D). Values 
represent mean ± SEM (error bars) from 
10-21 experiments. 
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Recovery from fast inactivation was measured at volt- 
ages ranging from —120 mV to — 70 mV (Fig. 5 A). At 
domain I SI, S151K slowed recovery and D152R was 
without effect (Fig. 5 B). Each of the remaining seven 
ENC charge-reversing mutations significantly acceler- 
ated recovery (Table 2). E161R (domain I S2), E598K 
(domain II SI), and N614K (domain II S2) produced 
an equivalent and moderate effect to accelerate recov- 
ery (Fig. 5 B). E1051R (domain III SI), D1069K (do- 
main III S2), E1373K (domain IV SI), and N1389K 
(domain TV S2) produced marked acceleration of channel 
recovery (Fig. 5, C and D). For domain III, these effects 
were charge dependent, as recovery time constants for 
E1051D (SI) and D1069E (S2) were similar to hNa v 1.4. 



In domain IV, E1373N and E1373D (SI) each accelerated 
recovery to a similar extent as for E1373K N1389E (S2) 
produced recovery kinetics similar to hNayl.4, and 
N1389D significandy slowed recovery at more hyperpolar- 
ized potentials, resulting in shallow voltage dependence. 

Deactivation kinetics were measured as delay in onset 
to recovery from fast inactivation (Fig. SI). Effects of 
ENC mutations on recovery delay paralleled those for 
recovery (Fig. S4 and Table 2) . For example, ENC muta- 
tion S151K (domain I SI), which slowed recovery, also 
prolonged recovery delay. ENC mutations that acceler- 
ated recovery shortened the delay to onset. Exceptions 
were mutations E161R and N1389D, which had no effect 
on recovery delay. Together, our findings indicate that 
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Figure 4. Fast inactivation kinet- 
ics in hNa v 1.4 and for ENC muta- 
tions. Time constants of current 
decay are shown for ENC muta- 
tions in domain I (A), domain 
II (B), domain III (SI, C; S2, D), 
and domain IV (SI, E; S2, F). Val- 
ues represent mean ± SEM (error 
bars) from 10-22 experiments. 
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domain III and domain IV ENC residues, as a conse- 
quence of their negative charge content, stabilize the 
fast-inactivated state. 

HCR mutations 

Putative countercharges in the bacterial VSM include 
SI residue N36 in NachBac (Yarov-Yarovoy et al., 2012) . 
The homologous residue N25 is observed in the crystal 
structure of NayAb as part of a hydrogen bond network 
with M21 in SI, S83 in S3, and R105 (R3) in S4 (Payandeh 
et al., 2011). We identified homologous residues in 
hNayl-4 (Fig. S3), and characterized the putative coun- 
tercharge in SI . To do this, we compared the effects of 
mutations at N144 (domain I), N591 (domain II), S1044 
(domain III) , and N1366 (domain IV) . These loci were 
considered to be the HCR. 

HCR mutations in domains I— III inhibited activation 
(Fig. 6 and Tables 1 and 2). Charge-reversing muta- 
tions N144R/K (domain I; Fig. 6, B and E) or N591R/K 



(domain II; Fig. 6, C and F) produced 16-17-mV depo- 
larizing shifts of Vo.5 and slowed channel activation. For 
domain I HCR mutations, effects were charge depen- 
dent, as N144D/E produced V0.5 values and activation 
kinetics similar to hNayl.4. In domain II, glutamate 
substitution at N591 rescued wild-type V0.5, whereas as- 
partate substitution produced a depolarizing shift of 
14 mV, similar to N591R/K. Both N591E and N591D 
partially restored wild-type activation kinetics. In do- 
main III, S1044R produced a depolarizing shift of 
12 mV, with S1044K/N producing a 5-mV shift (Fig. 6 D) . 
Activation was slowed by S1044R/K, but also by S1044D 
(Fig. 6 G). For the domain IV HCR locus N1366, aspar- 
tate substitution produced a significant, depolarizing 
shift of the conductance curve (Table 1). N1366D also 
decreased slope factor, whereas charge reversal (N1366R) 
increased slope factor. 

Figs. 7-9 show effects of HCR mutations on steady-state 
fast inactivation, kinetics of entry into, and recovery 
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Figure 5. Recovery kinetics in hNa v 1.4 and 
for ENC mutations. (A) Double pulse proto- 
col used to inactivate and recover channels. 
Traces are shown for recovery up to 50 ms, 
with some sweeps removed for clarity. Re- 
covery time constants are shown for ENC 
mutations in domain I and domain II (B), 
domain III (C), and domain IV (D). Values 
represent mean ± SEM (error bars) from 
10-18 experiments. 
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from that state. Lysine or arginine replacement at N144 
(domain I), N591 (domain II), and S1044 (domain III) 
produced a significant depolarizing shift of the h„ curve, 
as did arginine replacement at N1366 (domain IV; 
Fig. 7 and Table 1). N144D, S1044D, and N1336D/E 
each produced a significant hyperpolarizing shift of the 
h« curve. N591D and S1044N produced a depolarizing 
effect, and N144E and N591E produced h 0 .5 values simi- 
lar to hNayl.4. Slope factor was increased by N144D/E 
and by N591R, and decreased by N144K, N591K, and 
each of the mutations at N1366. 

Most charge-reversing mutations at the SI HCR locus 
slowed entry of channels into the fast-inactivated state 
(Fig. 8). N144R/K (domain I) and N591K (domain II) 
slowed fast inactivation, and both N144K and N591K 
reduced voltage dependence at voltages more positive 
than 0 mV. Charge substitutions at these loci produced 
inactivation kinetics similar to, or accelerated compared 
with, hNayl.4. In domains III and TV, charge-reversing 
mutations (S1044R, domain III; N1366R/K, domain IV) 



as well as charge-substituting mutations (S1044D, domain 
III; N1366D/E, domain IV) slowed entry of channels into 
the fast-inactivated state. N1366R slowed fast inactivation 
to the greatest extent across the voltage range tested. 
N1366D decreased the voltage dependence of fast inac- 
tivation, with accelerated kinetics at voltages more nega- 
tive than 0 mV. 

Recovery and its delay were most affected by HCR 
mutations in domains III and IV (Fig. 9, Fig. S5, and 
Table 2) . Addition of a positive charge at domain III 
HCR (S1044R) significantly accelerated recovery and 
its delay, whereas addition of a negative charge (S1044D) 
prolonged recovery and delay. For the domain IV HCR 
locus N1366, effects of mutations were less predictable. 
For example, N1366R slowed recovery, especially at 
hyperpolarized voltages, resulting in reduced voltage 
dependence. However, N1366R accelerated recovery 
delay, with increased voltage dependence. N1366E accel- 
erated recovery at depolarized voltages but prolonged 
recovery delay. 
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INC mutations 

In NachBac, residue E70, located in the inner portion 
of the S2 helix, sequentially interacts with positive charges 
R3 and R4 during activation (DeCaen et al., 2008, 
2009) . The homologous residue E59 in NayAb is ob- 
served as part of a hydrogen bond network between 
charged residues in S2, S3, and R108 (R4) in S4, as 
the INC (Payandeh et al., 2011). Alignments of bacte- 
rial sodium channel sequences with domain sequences 
from mammalian channels Nayl-2 (Yarov-Yarovoy et al., 
2012) or Na v 1.4 (Fig. 1; Zhang et al., 2012) show an- 
other conserved negative charge at the base of S3. Col- 
lectively, these INC residues are highly conserved in 
hNayl.4 (Fig. 1, A and B). Our alignment shows a gluta- 
mate one helical turn under the HCS region, in S2 (El 71, 
E624, E1079, E1399). Each of these negative charges in 
S2 is followed by a lysine (K175, K628, K1083, K1403) 
near the base of the helix. In S3, an inner residue, typi- 
cally asparagine (N194, N643, C1098, N1417) is fol- 
lowed by a conserved aspartate (D197, D646, D1101, 
D1420). We constructed charge-reversing mutations 
at each of these loci in S2 and S3, and characterized 
their effect on channel function. As neither D1101R nor 
D1101K expressed in oocytes, we tested valine replace- 
ment at this locus. 

Activation was inhibited by several mutations in the 
INC (Fig. 10 B and Tables 1 and 2). Overall, effects were 
notable only for charge reversals in the first two domains, 
and for residues closest to the hydrophobic region sur- 
rounding the gating pore. For example, charge reversal 
of INC residues at domain I S2 (E171R) or domain I S3 



(D197R) resulted in a depolarizing shift of V0.5 of 15 mV 
or more, but did not affect activation kinetics. Charge 
reversal at the homologous residues in domain II pro- 
duced significant depolarizing shifts of V0.5 (E624K, 
8 mV; D646K, 10 mV) and slowed activation rise time. 
Lesser depolarizing shifts in V0.5 were observed with mu- 
tation of the inner S3 residue in domain III (C1098K; 
5 mV) or in domain IV (N1417K; 6 mV); N1417K also 
slowed activation. Other mutations in the INC, which 
did not affect activation probability, did prolong rise 
time to peak activation of channels, including N194K 
(Dl S3), N643K (domain II S3), E1079R (domain III S2), 
E1399K (domain IV S2), and D1420K (domain IV S3). 

Several INC mutations affected one or both parame- 
ters of the hoo curve (Fig. 10C and Table 1). In domain I, 
E171R (S2) and D197R (S3) each produced a hyperpo- 
larizing shift of 10 mV, and N194K (S3) significandy 
reduced the slope factor, compared with wild-type chan- 
nels. In domain II, N643K (S3) produced a 5-mV hyper- 
polarizing shift of the h„ curve, and E624K (S2) increased 
the slope factor. In domain III, C1098K (S3) produced 
a 10-mV hyperpolarizing shift of h 0 .5 and decreased slope 
factor. Mutations at INC loci in domain IV decreased slope 
factor, and N14l7Kand D1420K (S3) produced a 4-mV 
hyperpolarizing shift. 

Charge-reversing INC mutations in each domain 
slowed entry of channels into the fast-inactivated state 
(Fig. 10, D and E; and Table 2). For domains I— III, ef- 
fects were largest for N194K (domain I S3) , D646K (do- 
main II S3), and E1079R (domain III S2), and N194K 
also reduced voltage dependence (Fig. 10 D). E171R 
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(domain I SI ) , D 197R (domain I S3), E624K and K628E 
(domain II S2), and N643K (domain II S3) slowed fast 
inactivation to a lesser extent. In domain IV, charge- 
reversing mutations at each of the four INC loci slowed 
fast inactivation (Fig. 10 E). E1399K (domain IV S2) and 
D1420K (domain IV S3), located near the gating pore 
constriction, slowed fast inactivation to a greater extent 
than K1403E and N1417K, located near the bases of S2 
and S3, respectively. 

INC mutations that produced 10 mV or greater hyper- 
polarizing shifts of the h*, curve (E171R, domain I S2; 
D197R, domain I S3; C1098K, domain III S3) slowed re- 
covery (Fig. 10, F and G; and Table 2). Each of these 
mutations also prolonged the delay in onset to recovery 
(Fig. S6 and Table 2). Several INC mutations that slowed 
the entry of channels into fast inactivation (N194K, do- 
main I S3; K1403E, domain IV S2; D1420K, domain IV S3) 
accelerated recovery, especially at depolarized voltages. 
D1420K slowed recovery at hyperpolarized voltages, re- 
sulting in decreased voltage dependence. 

DISCUSSION 

The gating particles proposed by Hodgkin and Huxley 
(1952) as the biophysical mechanism providing voltage 
sensitivity to excitable membranes have been identified 
as a repeating motif of positively charged residues in 
homologous S4 segments of the voltage-gated ion chan- 
nel superfamily (Yu et al., 2005). Results from experi- 
mental approaches using mutagenesis of S4 charge, 
site-specific toxins that trap S4 segments in resting or 
activated states, thiosulfonate reagents that discrimi- 
nate access of S4 residues to the intracellular or extra- 
cellular space, or fluorescence measurements of tagged 



S4 segments all support the premise that S4 segments 
move outward in response to membrane depolarization 
(for reviews see Bezanilla, 2008; Catterall, 2012). 

A challenge to sliding helix or helical screw models 
of S4 translocation is how voltage sensor movement is 
achieved in the interior of a highly hydrophobic mem- 
brane environment. Current models describe an S1-S4 
VSM that dictates the opening or closing of S5-S6 seg- 
ments comprising a pore module. This view of ion chan- 
nel gating is supported with experimental evidence in 
prokaryotic sodium channels that negative counter- 
charges in S1-S3 directly participate in S4 translocation 
(DeCaen et al., 2008, 2009, 201 1; Paldi and Gurevitz, 2010; 
Yarov-Yarovoy et al., 2012) , but has not yet been rigorously 
tested for eukaryotic sodium channels. 

The mechanisms underlying voltage gating in mam- 
malian sodium channels are crucial to our understand- 
ing of channelopathies in excitable tissues (George, 
2005; Cannon, 2006; Jurkat-Rott et al., 2010) . Investiga- 
tions of these disease-causing mutations have contrib- 
uted significandy to our understanding of the gating 
mechanisms underlying activation and fast inactivation. 
Crystal structures of bacterial channels have provided 
additional power toward the use of computational mod- 
els to suggest mechanisms of voltage gating (Vargas 
et al., 2012). Using the structural data provided by NayAb 
and NayRh, we created homology models of each do- 
main of the mammalian skeletal muscle sodium chan- 
nel and identified putative countercharges in SI, S2, 
and S3 segments oriented toward the S4 backbone. Our 
characterization of mutations at these loci show that 
activation is promoted by countercharges in domains 
I— III, and that sodium channel fast inactivation is dif- 
ferentially regulated by countercharges depending on 
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Figure 8. Kinetics of fast inactivation 
for HCR mutations. Plots show time 
constants for hNayl.4 and SI muta- 
tions in domain I (N144, A), domain II 
(N591, B), domain III (S1044, C), and 
domain IV (N1366, D). Values repre- 
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their location with respect to domain and the gating 
pore constriction. A view of domain-specific sodium 
channel functions is suggested in Fig. 11, showing the 
effects of VSM charge-reversing mutations. Activation 
is most affected by mutations in domains I and II, 
with fast inactivation most affected in domain IV. Do- 
main III mutations increase (INC) or decrease (HCR, 
ENC) steady-state fast inactivation. In general, muta- 
tions at ENC and HCR loci produce larger effects than 
those at INC loci. 

Negative countercharges facilitate activation of 
mammalian sodium channels 

In NachBac, transient pairings of residues E43 (SI) and 
D60 (S2) with S4 residues promote activation (DeCaen 
et al., 2008, 2009, 2011; Yarov-Yarovoy et al., 2012). We 
identified homologues in the VSM of each domain in 
hNayl .4. Our results strongly support a role in activation 



for SI and S2 ENC residues in domains I (S151, E161), 
II (E598, N614), and III (E1051, D1069). Charge-revers- 
ing mutations at these loci produced a depolarizing 
shift of Vo.s, and S151K, E161R, and E1051Ralso slowed 
activation kinetics. These results are consistent with the 
hypothesis that ENC residues in domains I— III interact 
with positively charged residues in S4 segments of these 
domains to stabilize S4 translocation in response to 
membrane depolarization during activation. 

The NachBac SI residue N36 has not yet been func- 
tionally characterized, but structural modeling suggests 
that this residue also participates as a countercharge 
in the VSM (Yarov-Yarovoy et al., 2012) . We identified 
homologous SI residues in hNayl. 4 located one heli- 
cal turn above the HCS (Fig. 1 and Fig. S3) . Our results 
show that hNayl-4 HCR SI residues contribute signi- 
ficantly to domain specificity of activation in hNayl-4. 
Mutational perturbation of activation was greatest for 
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vation for HCR mutations. (A) Traces 
showing inactivating and recovery 
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HCR residue N144 in domain I, and effects were charge 
dependent. Decreasing effects on activation were pro- 
duced by mutations of HCR residues N591, S1044, and 
N1366 in domains II— IV. The effect of HCR charge- 
reversing mutations in hNayl.4, inhibiting activation, 
are consistent with a model of the activation process 
in NacBac (Yarov-Yarovoy et al., 2012), showing that 
the homologous SI residue is part of a dynamic hydro- 



gen bond network that facilitates S4 translocation dur- 
ing activation. 

In the INC, activation probability was reduced by rever- 
sal of charge at El 71 (domain I S2),D197 (domainIS3), 
E624 (domain II S2), and D646 (domain II S3). Each 
of these mutations is located one helical turn under 
the HCS. Mutations of INC residues located two heli- 
cal turns under the HCS had no effect on activation. 
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Figure 1 0. Biophysical characterization of INC mutations. (A) Sodium currents in response to command depolarization to voltages 
from —90 mV to +60 mV. Effects of INC mutations are shown for activation (B), steady-state fast inactivation (C), entry into fast inacti- 
vation (D and E), and recovery (F and G). Legends identify residues and locus by domain and segment. Values represent mean ± SEM 
(error bars) from 10-23 experiments. 
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Collectively, our results indicate that ENC and HCR 
negative countercharges in domains I and II play the 
most critical role in activation, with contribution from 
domain III. 

Although the precise interactions of VSM counter- 
charges in hNayl-4 with S4 counterparts are uncertain, 
a comparison of our results to other studies suggests 
that the interaction of countercharge with outer positive 
charges in S4 is crucial to activation. Our free energy 
change calculations for the effect of point mutations of 
hNayl.4 countercharges yielded values of 1-3 kcal/mol, 
representing the energy of 1-2 hydrogen bonds in an 
a helix (Sheu et al., 2003). These values are similar to 
those for point mutations of NachBac S4 positive charges 
Rl, R3, and R4 (Yarov-Yarovoy et al., 2012), for which 



interaction of countercharge with S4 has been verified 
with disulfide locking experiments. Molecular dynamics 
simulations using homology models of rNayl.4 suggest 
interaction of Rl with a single ENC countercharge, and 
R2-R4 with ENC and INC countercharges during the ac- 
tivation process (Gosselin-Badaroudine et al., 2012). 
Given the relatively small change in activation free en- 
ergy associated with mutations of countercharge in 
hNayl .4, it is likely that interactions of S4 charge with 
residues nearest the gating pore constriction are the ac- 
tivation determinant. 

Each of the domain IV ENC mutations at SI (E1373) 
or at S2 (N1389) slowed activation without affecting con- 
ductance. This apparent conflict is most likely explained 
by the dramatic effect of these mutations to prolong 
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Figure 1 1 . Three-dimensional plot showing the effect of charge-reversing mutations of residues in ENC (green) , HCR (red) , and INC 
(blue) regions. Shown for each domain are shifts in inactivation kinetics at 20 mV (x axis), activation probability (y axis), and inactiva- 
tion probability (z axis). 
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fast inactivation. In response to depolarization, sodium 
channels activate and fast inactivate, with overlapping 
time courses at the macroscopic level (Aldrich et al., 
1983). Activation kinetics measured as time to peak acti- 
vation do not discriminate between slowed activation per 
se and prolonged time to peak current amplitude result- 
ing from impaired entry into fast inactivation. 

Simulations in rNayl.4 suggest interaction of counter- 
charges in S1-S3 with each of the first six positive charges 
in domain IV S4 during transitions between intermedi- 
ate states leading to activation (Gosselin-Badaroudine 
et al., 2012). In our homology models, these S4 residues 
are found in the ct helical region of S4, with K7 and R8 
located in the S4-S5 helix (unpublished data). From 
this computational data, and given the dramatic effect 
of domain IV ENC mutations on fast inactivation, it is 
likely that domain IV countercharges do interact with 
S4 residues during membrane depolarization, but that 
this interaction is not a determinant of pore opening. 
Our results indicating that negative countercharges in 
sodium channel domains I— III promote activation are 
consistent with previous studies showing that activation 
is dictated by S4 translocation in these domains (Chanda 
and Bezanilla, 2002; Armstrong, 2006) . 

Negative countercharges dictate fast inactivation gating 

Fast inactivation of the sodium channel and its subse- 
quent recovery determine the upper limit to action po- 
tential frequency in axon and muscle fiber membrane 
signaling. Our results show that negative countercharges 
are a crucial determinant in sodium channel fast inacti- 
vation and provide a greater understanding of its mech- 
anism. ENC residues in domains III (E1051, D1069) 
and IV (E1373, N1389) appear to play the most signifi- 
cant role in fast inactivation of hNayl.4. Charge reversal 
at each of these residues produced a depolarizing shift 
of the h„= curve, slowed entry of channels into fast in- 
activation, and accelerated their recovery. A previous 
study of domain IV S2 in the rat skeletal muscle sodium 
channel found that N1382C (Ma et al., 2009) produces 
the same pattern of effects on fast inactivation as we ob- 
served with charge reversal at the homologous locus in 
hNayl.4 (N1389K). 

Current models of sodium channel fast inactivation 
posit that voltage sensor translocation in response to 
membrane depolarization provides access for the in- 
activation particle to receptor sites in the S4-S5 linker 
of domain III (Smith and Goldin, 1997) and of domain 
IV (Lerche et al., 1997). Salt bridge formation between 
ENC residues in hNayl.4 and S4 residues during the 
depolarizing translocation preceding inactivation may 
explain the charge-dependent effects of domain III and 
domain IV ENC mutations on fast inactivation. How- 
ever, the molecular basis of hNayl.4 countercharge in- 
teraction with S4 remains to be shown experimentally. 
Additionally, our finding that E1373N (domain IV SI) 



reiterates each of the effects of E1373K on fast inactiva- 
tion indicates that the structure of the native residue 
is relevant. Length and orientation of the side chain 
of ENC countercharges glutamate, aspartate, and as- 
paragine are not equivalent and might confer specific 
constraints for salt bridge formation between these 
countercharges and putative S4 pairing partners in the 
native channel. 

Charge-reversing mutations in domain I produced 
moderate effects on parameters for fast inactivation. One 
interesting exception was S151K's effect of slowing recov- 
ery of channels from fast inactivation and to delay its 
onset. Assuming full absorption of S151K into the fast- 
inactivated state, effects on recovery are independent 
from slowed activation produced by this mutation. The 
prolonged delay in onset to recovery in S151K suggests 
that SI 51 promotes the deactivating transition of domain 
I S4, as channels leave the inactivated state and close 
(Kuo and Bean, 1994). A testable hypothesis is that this 
residue facilitates return of the nonimmobilized gating 
charge in that domain (Cha et al., 1999). 

In domain II, E598K and N614K produced effects 
on fast inactivation similar to those produced by muta- 
tions of ENC residues in domain III and domain IV, 
albeit lesser in magnitude. Voltage sensor mutations in 
all four domains of the sodium channel influence fast 
inactivation (Kontis and Goldin, 1997), and our results 
suggest that ENC residues in domain II, as well as in 
domain III and domain IV, play a role in fast inactiva- 
tion. However, it is not certain whether destabilization 
of the fast-inactivated state observed with E598K and 
N614K is a direct effect on fast inactivation per se. For 
example, domain II S4 gating transitions are coupled 
to other voltage sensor movements (Campos et al., 
2007b), and domain II ENC mutations might affect 
that coupling. 

HCR residues also appear to play a role in sodium 
channel fast inactivation. Charge reversal at N144, 
N591, SI 044, and N1366 decreased the probability that 
channels inactivate, and generally slowed the entry of 
channels into the fast-inactivated state. However, we 
found that predictions of charge substitution at these 
loci to produce wild-type fast inactivation parameters 
were not met, especially in terms of recovery and its de- 
lay. These results suggest that HCR residues may not 
influence fast inactivation by forming salt bridges with 
an S4-positive charge. An alternative possibility, hydro- 
gen bonding, is suggested from the crystal structure of 
NayAb (Payandeh et al. , 20 1 1 ) . In the closed-activated po- 
sition, S4 residue R105 (R3) makes hydrogen bonds with 
N25 (SI HCR), M29 (SI) , and S87 (S3). Homologous resi- 
dues in proximity to R3 are reiterated in the mammalian 
channel (Fig. S3). Further work is needed to investigate 
the specific targets of negative countercharges in the 
mammalian VSM, and structural modeling may help 
to clarify the roles of putative salt bridges or dynamic 
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hydrogen bonding in hNayl .4 gating transitions. None- 
theless, our results show that ENC and HCR residues 
are each important determinants of fast inactivation of 
the mammalian sodium channel. 

Several effects of INC mutations on fast inactivation are 
noteworthy. First, charge reversal or addition of positive 
charge at several loci produced a hyperpolarizing effect 
on the ho, curve, and none of the 12 mutations tested pro- 
duced the depolarizing shift that was consistendy observed 
with mutations of ENC or HCR loci. An interpretation of 
these effects is that at least some of the countercharges in 
the INC serve to restrict closed-state voltage sensor move- 
ment, whereas HCR and ENC charges facilitate S4 translo- 
cations, leading directly to pore opening. INC mutations 
that increased the likelihood of steady-state fast inactiva- 
tion also slowed channel recovery, suggesting that the INC 
contributes to stability of the fast-inactivated state through 
its impact on closed-state transitions. 

Each of four charge-reversing mutations in domain 
IV INC slowed entry of channels into the fast-inactivated 
state. A previous study showed that domain IV S2 muta- 
tion E1392C in rNa v 1.4 (Ma et al., 2009) slows fast inac- 
tivation, similar our finding with charge reversal at the 
homologous residue El 399. In a study of domain IV S3 
in hNayl. 4, cysteine substitution at D1420 slows fast in- 
activation to a greater extent than at N1417 (Nguyen and 
Horn, 2002), similar to our results with charge-reversing 
mutations D1420Kand N1417K. These findings support 
the hypothesis that intersegmental interactions of coun- 
tercharges with domain TV S4 are the crucial determinant 
of fast inactivation. However, our finding that domain I S3 
mutation N194K slows fast inactivation to the same ex- 
tent as observed for domain IV INC mutations suggests 
that countercharges in other domains also contribute 
to fast inactivation. 

Examination of crystal structures provided by bacterial 
sodium channels suggests that hydrogen bonding may be 
an important component of VSM interactions in the INC 
(Payandeh et al., 2011), in addition to the HCR. We ob- 
served a greater diversity of mutational perturbation of 
fast inactivation in each of these countercharge regions 
(HCR and INC) compared with the ENC. Future experi- 
mental work and structural modeling may provide ad- 
ditional insight into the mechanisms through which 
countercharges in each of these regions of the mamma- 
lian VSM (ENC, HCR, and INC) affect fast inactivation. 
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